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L'co/Lfir relations WITH CO ROTATIONAL LADDERS OF GALAXIES ACROSS THE HERSCHEL SPIRE 

ARCHIVE 

J. Kamenetzky^ N. Rangwala^, J. Glenn^, P. R. Maloney^, A. Conley'* 

ABSTRACT 

We present a catalog of all CO (J = 4—3 through J = 13—12), [Cl], and [Nil] lines available 
from extragalactic spectra from the Herschel SPIRE Fourier Transform Spectrometer (FTS) archive 
combined with observations of the low-J CO lines from the literature and from the Arizona Radio 
Observatory. This work examines the relationships between Lfir, Lqq, and Leo/Lco,i-o- We 
also present a new method for estimating probability distribution functions (PDFs) from marginal 
signal-to-noise ratio Herschel FTS spectra, which takes into account the instrumental “ringing” and 
the resulting highly correlated nature of the spectra. The slopes of log(LFiR) vs. log(LQQ) are 
linear for all mid- to high-J CO lines and slightly sublinear if restricted to (U)LlRGs. The mid- 
to high-J CO luminosity relative to CO J = 1-0 increases with increasing Lfir, indicating higher 
excitement of the molecular gas, though these ratios do not exceed ~ 180. For a given bin in Lfir, 
the luminosities relative to CO J=l—0 remain relatively flat from J = 6—5 through J= 13—12, across 
three orders of magnitude of Lfir. A single component theoretical photon-dominated region (PDR) 
model cannot match these flat SLED shapes, though combinations of PDR models with mechanical 
heating added qualitatively match the shapes, indicating the need for further comprehensive modeling 
of the excitation processes of warm molecular gas in nearby galaxies. 

Subject headings: galaxies: ISM - ISM: molecules - submillimeter: galaxies - submillimeter: ISM — 
surveys 


1. INTRODUCTION 

Within the multi-phase interstellar medium (ISM), 
molecular gas is the most intimately tied to star forma¬ 
tion, and therefore to the stellar lifecycle’s dramatic ef¬ 
fects on galaxy evolution. Though molecular hydrogen is 
the dominant component of such gas, pure H 2 rotational 
lines are difficult to detect and not particularly sensitive 
to the low temperatures of most molecular gas. Instead, 
(henceforth CO) and its isotopologues are used to 
trace the mass, kinematics, and excitation of molecular 
gas. The ground-level CO J = 1-0 line is widely used 
to estimate the total molecular mass in the interstellar 
medium, and ratios with higher lines provide informa¬ 
tion on the temperature and density of the emitting gas. 
Higher lines, however, are increasingly blocked by water 
absorption in Earth’s atmosphere. It was not until the 
launc h of the Herschel Space Observatory (Pilbratt et al. 
20101 that the CO ladder up to J = 13—12 was generally 
available for the ISM within our Galaxy and in nearby 
galaxies. 

Early SPIRE observations showed much brighter high- 
J CO emission than would be predicted by cool (T^in < 
50 K) molecular gas in giant molecular clouds, the type of 
gas respo nsible for the CO J = 1-0 and other low-J emis¬ 
sion (e.g., Panuzzo et al. 2010 Kanmne^kyet^^ 2012 


Rangwaia et al.||2011| Spinoglio et ^ 2012| Rigopoulou 

et al. 2013[ Perdra-Santaelia et al. 2013[ ). A warmer, 
denser (higher pressure) component ot molecular gas is 
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responsible for the emission of mid- (J = 4—3 to J = 
6—5) to high-J (J = 7—6 and above) CO lines (and even 
warmer emission can b e seen in much higher-J l i nes vi s- 
ible with PACS, as in Hailey-Dunsheath et al. (|2012 )). 
UV heating from young U and H stars creates Pliotodis- 
sociation Regions (PDRs), which can reproduce the ex¬ 
citation and emission of the low-J lines. However, PDR 
models often cannot explain the bright emission seen in 
high-J lines, which may require mechanical excitation 
via shocks, turbulence, winds, and other dynamical pro¬ 
cesses within galaxies. In addition to illuminating the 
excitation mechanisms of the gas, CO emission is also 


studied in the context of Kennicutt-Schmidt (Kennicutt 


1998 K-S law), which relate the gas surface density to 
the star formation rate (SFR) surface density. 

Now that Herschel's mission is complete, work is un¬ 
derway to e xamin e the full archival data set. [Kamenet^ 


zky et al. (2014) (henceforth K14) presented a two- 
component modeling procedure for a sample of galax¬ 
ies observed with the Herschel SPIRE Fourier Transform 
Spectrometer (FTS). In the 17 galaxy systems studied in 
that paper, the warm molecular gas accounted for about 
10% of the total molecular mass, but 90% of the CO 
luminosity. Here we expand this sample and compile a 
comprehensive, uniformly calibrated set of CO J = 4—3 
to J = 13-12, [Cl] (609 and 370 /rm), and [Nil] (205 
/rm) line fluxes for the galaxies observed by the Herschel 
SPIRE FTS, as well as a similarly matched set of CO 
J = l—0 to J = 3—2 lines from the literature and the 
Arizona Radio Observatory (ARO). A future paper will 
include a full, two-component likelihood analysis of each 
galaxy’s CO SLED, in order to derive cold and warm gas 
temperatures, densities, and masses, as in K14. 

The observations and processing are described in Sec¬ 
tion [^ Section [^ presents motivation for fitting the re- 
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lationships between L'qq and Lfir, our results broken 
down by subsamples, a nd comparison to th e si milar stud- 
ies of Lu et al. (2014), Greve et al. (2014) and Liu et al. 


(2015). Discussion of trends_and comparisons to theoret¬ 
ical niodels are in Sections IH 


2. OBSERVATIONS 


We compiled a list of successful extragalactic SPIRE 
FTS proposals (301 spectra) and searched the Herschel 
Science Archive (HSA) for the available data. In some 
cases, programs for higher or unknown redshift galaxies 
did not result in spectra with measurable CO emission, so 
those observations (about 74) are not presented. Table[^ 
lists the basic galaxy information and observation IDs for 
all galaxies for which at least one FTS line measurement 
or upper limit is reported]^ 


TABLE 1 

Galaxies and Observations Utilized 


Galaxy 

RA 

J2000 

Dec Log(LFiR) 

J2000 Lq 

Dl 

Mpc 

^det 

riul 

FTS Obsid 

Phot ObsID 

NGC0023 

0h09m53.36s 

+25d55m26.4s 

10.9 

68 

5 

7 

1342247622 

1342234681 

NGG34 

0hllm06.55s 

— 12d06m26.3s 

11.2 

85 

12 

0 

1342199253 

1342199383 

MCG-02-01-051 

0hl8m50.86s 

-10d22m37.5s 

11.2 

120 

8 

4 

1342247617 

1342234694 

IClO-Bll-1 

0h20m27.70s 

+59dl6m59.4s 

7.5 

1 

8 

5 

1342246982 

1342201446 

IRAS 00188-0856 

0h21m26.53s 

-08d39m27.1s 

12.2 

591 

4 

5 

1342246259 

1342234693 

ES0350-IG038=- 

0h36m52.46s 

-33d33ml7.4s 

10.8 

87 

1 

8 

1342246978 

1342199386 

NGC205-copeak 

0h40m24.10s 

+41d41m50.4s 

6.1 

1 

1 

9 

1342212315 

1342188661 

IRAS 00397-1312 

0h42ml5.53s 

-12d56m02.8s 

12.6 

1285 

0 

7 

1342246257 

1342234696 

NGG0232a 

0h42m45.82s 

-23d33m41.7s 

11.2 

95 

8 

4 

1342221707 

1342234699 

NGG253 

0h47m33.12s 

-25dl7ml7.6s 

10.3 

3 

12 

1 

1342210847 

1342199387 

I Zw !=■ 

0h53m34.94s 

+12d41m36.2s 

11.4 

272 

2 

9 

1342238246 

1342238252 

MCG-l-12-02-001 

0h54m03.61s 

+73d05mll.8s 

11.2 

72 

11 

2 

1342213377 

1342199365 

NGG0317B 

0h57m40.37s 

+43d47m32.4s 

11.0 

80 

8 

4 

1342239358 

1342238255 

IRAS 01003-2238 

Ih02m49.90s 

-22d21m57.3s 

11.9 

539 

1 

7 

1342246256 

1342234707 

3C 31 

Ih07m24.96s 

+32d24m45.2s 


75 

3 

5 

1342239344 

1342236245 

IC1623 

Ih07m47.00s 

-17d30m25.0s 

11.4 

86 

11 

1 

1342212314 

1342199388 

MCG-03-04-014=‘ 

Ihl0m08.92s 

-16d51mll.ls 

11.4 

152 

6 

6 

1342213442 

1342234709 

ESO244-G012‘" 

Ihl8m08.26s 

—44d27m43.0s 

11.1 

95 

5 

6 

1342221708 

1342234726 

GGCG436-030=>' 

Ih20m02.58s 

+14d21m42.5s 

11.5 

138 

9 

3 

1342213443 

1342237499 

ES0353-G020 

Ih34m51.29s 

-36d08ml5.0s 

10.8 

66 

7 

5 

1342247615 

1342234721 

IRASF01417-I-1651 

Ih44m30.52s 

+17d06m08.9s 


120 

7 

5 

1342239343 

1342237555 

NGG0695 

Ih51ml4.28s 

+22d34m55.2s 

11.4 

143 

6 

5 

1342224767 

1342238266 

Mrk 1014=- 

Ih59m50.21s 

+00d23m40.6s 

12.3 

763 

1 

7 

1342238707 

1342237540 

NGG0828 

2hl0m09.50s 

+39dllm24.7s 

11.1 

80 

4 

7 

1342239357 

1342239822 

NGG0877a 

2hl8m00.12s 

+14d32m34.2s 


57 

5 

7 

1342239342 

1342238267 

NGG 891-1 

2h22m33.41s 

+42d20m56.9s 

10.2 

10 

7 

5 

1342213376 

1342189430 

UGG01845 

2h24m07.89s 

+47d58mll.3s 

10.9 

70 

5 

7 

1342240022 

1342239799 

NGG0958 

2h30m42.80s 

-02d56m23.7s 

10.9 

82 

3 

7 

1342239339 

1342238277 

0235+164 

2h38m38.93s 

+16d36m59.3s 


4278 

0 

1 

1342249452 

1342224149 

NGG1068 

2h42m40.71s 

-00d00m47.8s 

10.9 

16 

13 

0 

1342213445 

1342189425 

NGG1056 

2h42m48.30s 

+28d34m27.1s 

9.7 

24 

5 

8 

1342204024 

1342226630 

UGG02238 

2h46ml7.50s 

+13d05m44.9s 

11.1 

93 

5 

5 

1342239340 

1342238270 

NGG1097 

2h46ml9.00s 

—30dl6m30.0s 

10.4 

16 

8 

5 

1342239337 

1342188586 

UGG02369 

2h54m01.81s 

+14d58ml4.3s 

11.4 

142 

8 

3 

1342239341 

1342239831 

NGG1222 

3h08m56.74s 

-02d57ml8.5s 

10.4 

35 

8 

5 

1342239354 

1342239262 

UGG02608 

3hl5m01.24s 

+42d02m09.2s 

11.1 

104 

4 

5 

1342239356 

1342239819 

NGG1266 

3hl6m00.70s 

-02d25m38.0s 

10.2 

31 

13 

0 

1342239353 

1342189424 

IRAS 03158+4227“^ 

3hl9ml2.40s 

+42d38m28.0s 

12.4 

623 

4 

4 

1342224764 

1342226656 

3C 84 

3hl9m48.16s 

+41d30m42.1s 

10.8 

78 

11 

2 

1342249054 

1342203614 

NGG1365-SW 

3h33m35.90s 

-36d08m35.0s 

10.8 

21 

10 

2 

1342204021 

1342201432 

NGG1365-NE 

3h33m36.60s 

—36d08m20.0s 

10.8 

21 

10 

3 

1342204020 

1342201432 

NGG1377 

3h36m39.10s 

-20d54m08.0s 

9.7 

24 

12 

1 

1342239352 

1342189505 

NGG1482 

3h54m38.90s 

-20d30m09.0s 

10.5 

25 

12 

1 

1342248233 

1342189504 

IRAS 03521+0028“^ 

3h54m42.19s 

+00d37m02.0s 

12.3 

709 

2 

7 

1342238704 

1342239850 

UGG02982 

4hl2m22.53s 

+05d32m50.4s 

10.9 

77 

3 

7 

1342240021 

1342239938 

ES0420-G013“" 

4hl3m49.65s 

-32d00m24.1s 

10.7 

49 

8 

5 

1342242590 

1342227719 

NGG1572“‘ 

4h22m42.81s 

—40d36m03.2s 

11.0 

86 

6 

5 

1342242588 

1342227720 

IRAS04271+3849 

4h30m33.10s 

+38d55m48.4s 

10.9 

86 

6 

4 

1342227786 

1342229106 

NGG1614 

4h33m59.85s 

—08d34m44.0s 

11.3 

68 

12 

1 

1342192831 

1342203628 

UGG03094 

4h35m33.75s 

+19dl0ml7.5s 

11.1 

108 

4 

6 

1342227522 

1342239944 

MCG-05-12-006“" 

4h52m04.96s 

-32d59m25.9s 

10.9 

78 

6 

4 

1342242589 

1342229237 

IRAS F05189-2524“" 

5h21m01.47s 

-25d21m45.4s 

11.8 

185 

10 

2 

1342192833 

1342203632 

IRAS05223+1908 

5h25ml6.65s 

+19dl0m48.5s 


130 

2 

0 

1342228738 

1342229652 

MCG+08-11-002 

5h40m43.65s 

+49d41m41.8s 

11.2 

86 

9 

3 

1342230414 

1342229112 

NGG1961 

5h42m04.37s 

+69d22m41.9s 

10.7 

61 

9 

4 

1342228708 

1342227742 

UGG03351 

5h45m48.00s 

+58d42m03.7s 

11.1 

67 

6 

6 

1342230415 

1342229115 

IRAS05442+1732 

5h47mll.l5s 

+17d33m47.2s 

11.0 

81 

5 

5 

1342230413 

1342229653 


The spectra for this survey came from the followingprograms, 
with the total number of observations presented in Tableluin paren¬ 
theses: DT1 ji1u_ 1 (92), KPOT_pvanderw_l (31), 0Tl_dfarrah_l 
(28), 0Tl_jsinith01_l (23), 0Tl_pogle01_l (13), KPGT_cwilso01_l 


(12), 0Tl_lyoung_l (8), GTl_lspinogl_2 (10), 0Tl_pvaiiderw_4 (5) 
GT2_vleboute_3 (2), QT2_vkulkarn_3 (2), KPGT_rguesten_l (1) 
0Tl_dmarrone_l (0), 0Tl_rivlson_l (0), 0T2_drigopou_3 (0) 
0T2_rivison_2 (0). 
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Lqo/^fir Relations Across the Herschel SPIRE Archive 

TABLE 1 — Continued 


Galaxy 

RA 

J2000 

Dec 

J2000 

Log(LFia) 

L© 

Dl 

Mpc 

'^det 


FTS Obsid 

Phot ObsID 

IRAS 06035-7102 

6h02m54.01s 

-71d03ml0.2s 

12.0 

353 

10 

2 

1342230420 


UGC03410a 

6hl4m29.64s 

+80d26m59.4s 

10.8 

61 

3 

8 

1342231072 

1342229131 

NGC2146-NW 

6hl8m36.70s 

+78d21m32.0s 

10.8 

17 

12 

0 

1342219554 

1342191186 

NGC2146-nuc 

6hl8m38.60s 

+78d21m24.0s 

10.8 

17 

11 

2 

1342204025 

1342191186 

NGC2146-SE 

6hl8m40.50s 

+78d21ml6.0s 

10.8 

17 

11 

2 

1342219555 

1342191186 

IRAS 06206-6315=- 

6h21m01.21s 

-63dl7m23.5s 

12.0 

411 

4 

5 

1342231038 

1342226638 

ES0255-IG007= 

6h27m21.63s 

-47dl0m36.3s 


166 

7 

4 

1342231084 

1342226643 

UGC03608 

6h57m34.42s 

+46d24ml0.7s 

11.1 

97 

6 

6 

1342228744 

1342229649 

NGC2342b 

7h09ml2.09s 

+20d36ml3.1s 

10.8 

77 

5 

6 

1342228730 

1342230778 

NGC2342a 

7h09ml8.05s 

+20d38ml0.0s 

10.8 

77 

5 

5 

1342228729 

1342230778 

NGG2369 

7hl6m37.60s 

-62d20m35.9s 

10.8 

43 

10 

3 

1342231083 

1342229670 

NGC2388a 

7h28m53.43s 

+33d49m08.4s 

11.0 

62 

7 

6 

1342231071 

1342229477 

MGG-l-02-20-003 

7h35m43.44s 

+Ild42m34.8s 

10.8 

72 

7 

6 

1342228728 

1342229463 

IRAS 07598+6508= 

8h04m30.45s 

+64d59m52.2s 

12.1 

693 

0 

10 

1342253659 

1342229642 

B2 0827+24 

8h30m52.09s 

+24dl0m59.8s 


5818 

0 

1 

1342253660 

1342230773 

IRAS 08311-2459= 

8h33m20.60s 

-25d09m33.7s 

12.2 

451 

9 

1 

1342230421 

1342230796 

He2-10 

8h36ml5.18s 

—26d24m33.9s 

9.6 

10 

9 

3 

1342245083 

1342196888 

IRAS08355-4944 

8h37m01.86s 

-49d54m30.0s 


110 

7 

5 

1342231975 

1342226978 

NGC2623 

8h38m24.08s 

+25d45ml6.6s 

11.4 

81 

12 

0 

1342219553 

1342206174 

IRAS 08572+3915 

9h00m25.39s 

+39d03m54.4s 

11.8 

261 

2 

8 

1342231978 

1342230749 

IRAS09022-3615 

9h04ml2.72s 

-36d27m01.3s 

12.0 

262 

11 

1 

1342231063 

1342230799 

NGC2764 

9h08ml7.47s 

+21d26m36.0s 

10.0 

40 

5 

7 

1342231057 

1342245567 

NGC2798 

9hl7m22.90s 

+41d59m59.0s 

10.4 

28 

12 

1 

1342252892 

1342197287 

UGC05101 

9h35m51.65s 

+61d21mll.3s 

11.8 

176 

9 

3 

1342209278 

1342204962 

NGC2976_00 

9h47m07.84s 

+67d55m52.3s 


4 

1 

10 

1342228706 

1342192106 

M81 

9h55m33.17s 

+69d03m55.0s 

9.2 

4 

3 

9 

1342209851 

1342185538 

M82 

9h55m52.22s 

+69d40m46.9s 

10.4 

4 

13 

0 

1342208389 

1342185537 

NGC3077 

10h03ml9.10s 

+68d44m02.0s 

7.7 

1 

1 

9 

1342228745 

1342193015 

NGC3110a 

10h04m02.09s 

—06d28m28.6s 

11.0 

73 

4 

6 

1342231971 

1342234843 

3C 236 

10h06m01.74s 

+34d54ml0.4s 


451 

0 

6 

1342246988 

1342246613 

NGC3221 

10h22m20.20s 

+21d34m22.4s 

10.7 

61 

2 

8 

1342221714 

1342246610 

NGC3227= 

10h23m30.58s 

+19d51m54.2s 

9.7 

18 

12 

1 

1342209281 

1342197318 

NGC3256 

10h27m51.27s 

-43d54ml3.8s 

11.3 

38 

13 

0 

1342201201 

1342200126 

IRAS 10378+1109 

10h40m29.17s 

+10d53ml8.3s 

12.1 

631 

2 

8 

1342247118 

1342234867 

ESO264-G036 

10h43m07.68s 

—46dl2m44.9s 

10.9 

89 

6 

4 

1342249044 

1342236204 

NGC3351 

10h43m57.70s 

+Ild42ml4.0s 

9.7 

13 

9 

4 

1342247117 

1342198885 

ESO264-G057 

10h59m01.82s 

-43d26m25.9s 

10.7 

72 

6 

5 

1342249043 

1342236203 

IRASF10565+2448 

10h59ml8.17s 

+24d32m34.4s 

11.8 

192 

10 

2 

1342247096 

1342234869 

NGG3521 

Ilh05m48.60s 

—00d02m09.0s 

10.1 

12 

6 

3 

1342247743 

1342198568 

IRAS 11095-0238 

Ilhl2m03.38s 

-02d54m23.8s 

12.0 

482 

4 

5 

1342247760 

1342234863 

NGC3627 

Ilh20ml5.00s 

+12d59m30.0s 

10.2 

12 

13 

0 

1342247604 

1342198883 

NGC3665 

Ilh24m43.67s 

+38d45m46.0s 

9.7 

32 

2 

7 

1342247121 

1342222667 

Arp299-B 

Ilh28m31.00s 

+58d33m41.0s 

11.6 

49 

13 

0 

1342199249 

1342199345 

Arp299-G 

Ilh28m31.00s 

+58d33m50.0s 

11.6 

49 

13 

0 

1342199250 

1342199345 

Arp299-A 

Ilh28m33.63s 

+58d33m47.0s 

11.6 

49 

13 

0 

1342199248 

1342199345 

PG 1126-041 

Ilh29ml6.66s 

-04d24m07.6s 


266 

1 

9 

1342247119 

1342247271 

ESO 320-G030 

llh53mll.72s 

-39d07m48.9s 

11.0 

45 

12 

1 

1342210861 

1342200129 

NGC3982= 

Ilh56m28.13s 

+55d07m30.9s 

9.8 

21 

5 

5 

1342209277 

1342186862 

NGG4038 

12h01m53.00s 

-18d52m01.0s 


23 

6 

5 

1342210860 

1342188686 

N GG4038overlap 

12h01m54.90s 

-18d52m46.0s 


23 

8 

4 

1342210859 

1342188686 

NGG4051= 

12h03m09.61s 

+44d31m52.8s 

9.5 

14 

9 

4 

1342209276 

1342210502 

IRAS 12071-0444 

12h09m45.12s 

-05d01ml3.9s 

12.1 

591 

2 

8 

1342248239 

1342234858 

NGG4151= 

12hl0m32.58s 

+39d24m20.6s 


18 

5 

7 

1342209852 

1342188588 

NGG4194 

12hl4m09.63s 

+54d31m36.1s 

10.7 

39 

10 

3 

1342231069 

1342230869 

IRAS12116-5615 

12hl4m22.17s 

-56d32m32.8s 

11.3 

115 

10 

2 

1342249462 

1342226974 

NGG4254 

12hl8m49.60s 

+14d24m59.0s 

10.8 

36 

5 

4 

1342236997 

1342187173 

NGG4321 

12h22m54.90s 

+15d49m21.0s 

10.4 

25 

8 

4 

1342247572 

1342187322 

NGG4388 

12h25m46.75s 

+12d39m43.5s 

10.4 

38 

10 

3 

1342210849 

1342248482 

NGG4459 

12h29m00.03s 

+13d58m42.8s 

9.1 

19 

3 

7 

1342248411 

1342200118 

NGG4526 

12h34m03.03s 

+07d41m56.9s 

9.1 

10 

6 

6 

1342224762 

1342234889 

NGG4536 

12h34m27.00s 

+02dllml7.0s 

10.5 

27 

9 

4 

1342237025 

1342189455 

NGG4569 

12h36m49.80s 

+13d09m46.0s 

7.5 

1 

10 

3 

1342248251 

1342188777 

TOL1238-364 

12h40m52.85s 

-36d45m21.1s 

10.4 

46 

7 

6 

1342213381 

1342202200 

NGG4631 

12h42m08.00s 

+32d32m29.0s 

10.4 

12 

7 

3 

1342247573 

1342188756 

NGG4710 

12h49m38.96s 

+15d09m55.8s 

9.6 

18 

5 

8 

1342247120 

1342188766 

NGG4736 

12h50m53.00s 

+41d07ml4.0s 

9.9 

8 

7 

4 

1342245851 

1342188754 

Mrk 231 

12h56ml4.23s 

+56d52m25.2s 

12.2 

188 

11 

1 

1342210493 

1342201218 

NGG4826 

12h56m43.70s 

+21d40m58.0s 

9.7 

9 

8 

3 

1342246992 

1342188764 

MGG-02-33-098 

13h02ml9.80s 

-15d46m03.5s 

10.7 

69 

4 

7 

1342247567 

1342234810 

ES0507-G070 

13h02m52.34s 

-23d55ml7.8s 

11.2 

91 

11 

1 

1342248421 

1342234813 

NGC5010 

13hl2m26.39s 

-15d47m51.7s 

10.6 

43 

6 

5 

1342236996 

1342234809 

IRAS13120-5453 

13hl5m06.35s 

-55d09m22.7s 

12.0 

132 

12 

0 

1342212342 

1342226970 

NGC5055 

13hl5m49.30s 

+42d01m45.0s 

10.1 

11 

5 

6 

1342237026 

1342188753 

Arpl93 

13h20m35.34s 

+34d08m22.2s 

11.4 

105 

12 

0 

1342209853 

1342198191 

NGG5104 

13h21m23.09s 

+00d20m33.3s 

10.9 

82 

5 

6 

1342247566 

1342236168 

MGG-03-34-064= 

13h22m24.42s 

-16d43m42.7s 


74 

8 

3 

1342249041 

1342236178 
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Galaxy 

RA 

J2000 

Dec 

J2000 

Log(LFia) 

L© 

Dl 

Mpc 

'^det 


FTS Obsid 

Phot ObsID 

Cen A 

13h25m27.61s 

-43d01m08.8s 

9.7 

4 

8 

5 

1342204037 

1342188663 

NGC5135 

13h25m44.06s 

-29d50m01.2s 

11.0 

58 

12 

0 

1342212344 

1342202248 

ESO 173-G015 

13h27m23.78s 

-57d29m22.2s 

11.2 

39 

13 

0 

1342202268 

1342203562 

NGG5194 

13h29m52.71s 

+47dllm42.6s 

10.4 

11 

7 

6 

1342201202 

1342188589 

IC4280 

13h32m53.35s 

-24dl2m25.4s 

10.7 

70 

6 

4 

1342249042 

1342236191 

M83 

13h37m00.92s 

-29d51m56.7s 

10.4 

7 

10 

2 

1342212345 

1342188664 

Mrk 273 

13h44m42.11s 

+55d53ml2.7s 

12.0 

168 

11 

1 

1342209850 

1342201217 

4C 12.50=" 

13h47m33.36s 

+12dl7m24.2s 

12.0 

561 

1 

9 

1342237024 

1342234792 

UGC08739 

13h49ml4.28s 

+35dl5ml9.8s 

10.8 

76 

6 

6 

1342247123 

1342236144 

ES0221-IG010 

13h50m56.87s 

-49d03ml8.5s 

10.5 

39 

5 

6 

1342249461 

1342238293 

Mrk 463“* 

13h56m02.87s 

+18d22ml9.5s 

11.2 

226 

4 

7 

1342249047 

1342236151 

M101_02 

14h03m41.36s 

+54dl9m04.9s 

10.1 

8 

2 

10 

1342230417 

1342188750 

OQ 208“* 

14h07m00.39s 

+28d27ml4.7s 


348 

3 

6 

1342247769 

1342234785 

NGG5653 

14h30m09.88s 

+31dl2m56.3s 

10.8 

55 

5 

8 

1342247565 

1342236146 

IRAS 14348-1447 

14h37m38.26s 

-15d00m24.6s 

12.1 

371 

8 

3 

1342249457 

1342238301 

NGG5713 

14h40ml 1.50s 

—00dl7m20.0s 

10.5 

29 

8 

5 

1342248250 

1342189520 

IRAS 14378-3651 

14h40m59.01s 

-37d04m32.0s 

11.9 

303 

10 

1 

1342227456 

1342238295 

Mrk 478“* 

14h42m07.46s 

+35d26m22.9s 

11.1 

358 

0 

9 

1342238710 

1342238333 

NGG5734a 

14h45m08.98s 

-20d52ml3.4s 

10.7 

59 

3 

10 

1342248417 

1342227731 

3C 305“* 

14h49m21.80s 

+63dl6ml5.3s 


187 

1 

4 

1342236998 

1342234915 

VV340a“* 

14h57m00.66s 

+24d37m05.1s 


145 

6 

5 

1342238241 

1342234779 

IC4518ABa 

14h57m41.15s 

-43d07m56.2s 


68 

7 

4 

1342250514 

1342239895 

NGG5866“* 

15h06m29.50s 

+55d45m47.6s 

9.4 

14 

7 

4 

1342238708 

1342188749 

CGCG049-057 

15hl3ml3.09s 

+07dl3m31.8s 

11.1 

59 

11 

2 

1342212346 

1342203077 

3C 315 

15hl3m40.08s 

+26d07m31.2s 


498 

0 

4 

1342239350 

1342234777 

VV705“* 

15hl8m06.13s 

+42d44m44.5s 


181 

10 

2 

1342238712 

1342229532 

ES0099-G004 

15h24m57.99s 

-63d07m30.2s 

11.4 

125 

10 

2 

1342230419 

1342229209 

IRAS 15250+3609 

15h26m59.40s 

+35d58m37.5s 

11.8 

248 

2 

9 

1342238711 

1342234775 

NGG5936 

15h30m00.80s 

+12d59m21.7s 

10.8 

61 

6 

5 

1342249046 

1342238324 

Arp220 

15h34m57.12s 

+23d30mll.5s 

12.0 

81 

11 

2 

1342190674 

1342188687 

NGC5990 

15h46ml6.40s 

+02d24m54.7s 

10.7 

57 

6 

7 

1342240016 

1342238312 

IRAS 15462-0450 

15h48m56.81s 

-04d59m33.6s 

12.0 

456 

3 

7 

1342249045 

1342238307 

3C 326 

15h52m09.07s 

+20d05m48.4s 


407 

0 

7 

1342250516 

1342238327 

PKS 1549-79 

15h56m58.87s 

-79dl4m04.3s 


690 

0 

8 

1342253671 

1342239890 

NGC6052 

16h05ml2.94s 

+20d32m36.9s 

10.8 

71 

6 

6 

1342212347 

1342229560 

IRAS 16090-0139 

16hllm40.48s 

-01d47m05.6s 

12.3 

618 

6 

4 

1342238699 

1342229565 

PG 1613+658 

16hl3m57.18s 

+65d43m09.6s 

11.5 

600 

0 

8 

1342242593 

1342238336 

CGGG052-037 

16h30m56.60s 

+04d04m58.3s 

11.1 

109 

8 

2 

1342251284 

1342229572 

NGC6156 

16h34m52.50s 

-60d37m07.7s 

10.8 

45 

8 

4 

1342231041 

1342229213 

ES0069-IG006 

16h38ml 1.84s 

-68d26m08.5s 

11.7 

203 

7 

4 

1342231040 

1342230810 

IRASF16399-0937 

16h42m40.10s 

—09d43ml3.6s 

11.3 

118 

8 

4 

1342251334 

1342229188 

NGC6240 

16h52m58.89s 

+02d24m03.4s 

11.6 

108 

13 

0 

1342214831 

1342203586 

IRASF16516-0948 

16h54m23.81s 

-09d53m21.4s 

11.0 

100 

6 

5 

1342251335 

1342229189 

NGC6286b 

16h58m23.99s 

+58d57m20.3s 

11.1 

85 

2 

8 

1342231068 

1342229148 

NGC6286a 

16h58m31.56s 

+58d56ml2.2s 

11.1 

85 

7 

3 

1342221715 

1342229148 

IRASF17138-1017 

17hl6m35.82s 

-10d20m41.5s 

11.1 

76 

8 

3 

1342230418 

1342229190 

IRAS F17207-0014 

17h23m21.96s 

-00dl7m00.9s 

12.2 

190 

11 

1 

1342192829 

1342203587 

ESO138-G027 

17h26m43.30s 

-59d55m55.6s 

11.1 

88 

6 

6 

1342231042 

1342229216 

UGC11041 

17h54m51.82s 

+34d46m34.3s 

10.8 

74 

4 

5 

1342231061 

1342229169 

IRAS17578-0400 

18h00m31.86s 

-04d00m53.3s 

11.2 

62 

8 

4 

1342231047 

1342229187 

NGC6621 

18hl2m55.31s 

+68d21m46.8s 

11.0 

92 

5 

4 

1342221716 

1342220865 

IC4687 

18hl3m39.63s 

-57d43m31.3s 

11.1 

73 

10 

3 

1342192993 

1342204955 

IRAS F18293-3413 

18h32m41.13s 

-34dllm27.5s 

11.5 

78 

11 

2 

1342192830 

1342204954 

IC4734 

18h38m25.60s 

-57d29m25.1s 

11.0 

67 

8 

4 

1342240013 

1342229222 

NGC6701 

18h43ml2.56s 

+60d39mll.3s 

10.9 

62 

9 

4 

1342231994 

1342229137 

IRAS 19254-7245“* 

19h31m20.50s 

-72d39m21.8s 

11.8 

270 

7 

2 

1342231039 

1342206210 

IRAS 19297-0406 

19h32m22.00s 

—04d00m02.0s 

12.2 

387 

5 

7 

1342231078 

1342230837 

ESO339-G011 

19h57m37.59s 

-37d56m08.5s 

10.8 

82 

3 

8 

1342231990 

1342230821 

3C 405 

19h59m28.36s 

+40d44m01.9s 


252 

2 

7 

1342246994 

1342230853 

IRAS 20087-0308 

20hllm23.87s 

-02d59m50.7s 

12.2 

480 

7 

3 

1342231049 

1342230838 

IRAS 20100-4156 

20hl3m29.54s 

—41d47m34.9s 

12.4 

595 

6 

4 

1342245106 

1342230817 

MGG+04-48-002a 

20h28m35.02s 

+25d44m00.6s 

10.9 

65 

4 

6 

1342221682 

1342233320 

NGC6926 

20h33m06.08s 

-02d01m38.7s 

11.0 

87 

4 

7 

1342231050 

1342218992 

NGC6946 

20h34m52.30s 

+60d09ml4.0s 

9.8 

5 

11 

1 

1342243603 

1342188786 

NGC6946_05 

20h35ml2.01s 

+60d08m55.2s 

9.8 

5 

4 

8 

1342224769 

1342188786 

IRAS 20414-1651 

20h44ml8.21s 

— 16d40ml6.2s 

12.0 

392 

3 

8 

1342243623 

1342231345 

3C 424 

20h48ml2.03s 

+07d01ml7.5s 


586 

0 

9 

1342255797 

1342244149 

IC 5063“* 

20h52m02.10s 

-57d04m06.6s 

10.2 

46 

3 

8 

1342242619 

1342206208 

CGGG448-020“* 

20h57m24.33s 

+17d07m38.3s 

11.7 

161 

10 

2 

1342221679 

1342233327 

ESO286-IG019 

20h58m26.79s 

—42d39m00.6s 

11.8 

185 

11 

1 

1342245107 

1342230815 

ESO286-G035 

21h04mll.l3s 

-43d35m34.1s 

10.8 

73 

7 

5 

1342216901 

1342230813 

3C 433 

21h23m44.60s 

+25d04m27.1s 


465 

0 

9 

1342245864 

1342234675 

NGG7130 

21h48ml9.50s 

-34d57m04.7s 

11.1 

69 

12 

1 

1342219565 

1342210527 

NGG7172 

22h02m01.91s 

-31d52mll.3s 

10.2 

36 

4 

7 

1342219549 

1342209301 

ESO467-G027 

22hl4m39.85s 

-27d27m50.5s 

10.8 

74 

2 

8 

1342245108 

1342245428 

IC5179 

22hl6m09.13s 

-36d50m36.6s 

10.9 

48 

6 

4 

1342245109 

1342244158 
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Galaxy 

RA 

J2000 

Dec 

J2000 

Log(LFiR) 

L© 

Dl 

Mpe 

^det 


ETS Obsid 

Phot ObsID 

NGC7331 

22h37m04.10s 

+34d24m56.0s 

10.3 

15 

6 

4 

1342245871 

1342189532 

UGG12150 

22h41ml2.19s 

+34dl4m56.2s 

11.1 

96 

6 

5 

1342221699 

1342220870 

IRAS 22491-1808 

22h51m49.26s 

-17d52m23.5s 

11.9 

345 

9 

2 

1342245082 

1342234671 

NGG7465 

23h02m00.96s 

+15d57m53.4s 

9.7 

30 

3 

8 

1342245869 

1342234763 

NGG7469 

23h03ml5.62s 

+08d52m26.4s 

11.3 

72 

13 

0 

1342199252 

1342196915 

ES0148-IG002 

23hl5m46.72s 

-59d03ml5.1s 

11.8 

193 

11 

1 

1342245110 

1342209299 

IC5298 

23hl6m00.64s 

+25d33m23.7s 

11.3 

122 

7 

5 

1342221700 

1342234766 

NGG7552 

23hl6ml0.77s 

-42d35m05.4s 

10.7 

21 

12 

1 

1342198428 

1342210528 

NGG7591 

23hl8ml6.26s 

+06d35m08.8s 

10.8 

72 

6 

7 

1342257346 

1342234758 

NGG7592^ 

23hl8m22.08s 

-04d24m57.6s 

11.1 

106 

7 

4 

1342221702 

1342234750 

NGG7582 

23hl8m23.50s 

-42d22ml4.0s 

10.5 

21 

11 

2 

1342209280 

1342210529 

IRAS 23230-6926 

23h26m03.62s 

-69dl0ml8.8s 

12.1 

482 

6 

4 

1342246276 

1342230806 

NGG7674 

23h27m56.68s 

+08d46m43.6s 

11.2 

130 

5 

7 

1342245858 

1342234929 

IRAS 23253-5415 

23h28m06.10s 

-53d58m31.0s 

12.1 

595 

4 

6 

1342246277 

1342234737 

NGG7679a^ 

23h28m46.61s 

+03d30m41.8s 

10.8 

75 

6 

7 

1342221701 

1342234755 

IRAS 23365+3604 

23h39m01.27s 

+36d21m08.7s 

12.0 

290 

7 

5 

1342224768 

1342234919 

NGG7771 

23h51m24.88s 

+20d06m42.6s 

11.1 

63 

10 

3 

1342212317 

1342199379 

Mrk331 

23h51m26.80s 

+20d35m09.9s 

11.2 

81 

12 

0 

1342212316 

1342234682g 


Note. — firmdet and riui indicate the number of 3 sigma de¬ 
tections and upper limits, respectively, reported in Table 1^ out of 
our 13 fitted lines: CO J=4—3 to J= 13—12, two [Cl] lines, and 
one [Nil], 

^ Indicates photometry correction was not performed on extended 
FTS spectrum, see Section|2.2| 


2 . 1 . Herschel SPIRE Photometry and FTS Spectra 

All spectra used in the sample were repro¬ 
cessed with HlPhjj developer’s version 13.0.3849 and 
spire_cal_13_0, ootained from Rosalind Hopwood on 
2014 September 30. This calibration corrects for rapidly 
changing telescope temperatures near the beginning of 
observat ion cycles, which has the largest effect on faint 
sources (Swinyard et al. 20141. Overall the calibration 
errors, even from earlier calibration sets, are within 6 % 
for point sources and 7% for extended sources. The ma¬ 
jority of the observations were done in sparse sampling 
mode, for which we took the spectra from the central pix¬ 
els (SLWC3, SSWD4). For those mapping observations 
in intermediate or fully sampled mode, we extracted the 
spectrum from the pixel correspondiM to the central co¬ 
ordinates of the map (those in Table IR. All SPIRE pho¬ 
tometry was downloaded in 2014 September from the 
Herschel Science Archive (SPG vll.1.0) and not repro¬ 
cessed. 


2 . 2 . Source/Beam Correction 

The majority of the sources in our band are quite point¬ 
like compared to the SPIRE FTS beam, which varies 
from ~ 45" to 17". Because the beam size is discontin¬ 
uous between the upper frequency end of the SLW band 
and the lower frequency end of the SSW band, galaxies 
which are not point-like will show a notable discontinuity 
at this juncture. Even galaxies that are relatively small 
compared to the beam, but still not perfectly point-like, 
will show this discontinuity and require a correction to 
properly compare the emission across the SPIRE band¬ 
pass. This is necessary because we only use the central 
FTS detectors. An example is shown in Appendix [A| 

® Hess, HSpot, and HIPE are joint developments by the Her¬ 
schel Science Ground Segment Consortium, consisting of ESA, the 
NASA Herschel Science Center, and the HIFI, PACS and SPIRE 
consortia. 


We perform the same source/beam correction as de¬ 
scribed in full detail in Section 2.2 of K14. Briefly, we 
use the SPIRE Photometer Short Wave (PSW, 250 /tm) 
maps, observed with 19732 beams, and convolve them to 
larger beam sizes (B{,) and measure the new peak flux 
density. We compare this flux density to that of a 6 = 
4375 beam, which corresponds to the beam at the CO J = 
4—3 transition. The ratio of the two flux densities, ^ 76 , 43 . 5 , 
as a function of beam size is between that expected for a 
point source ( 1 ) and a fully extended source (r2t,/B43 5). 
For each galaxy’s unique distribution, for any beam size, 
we have a value of ? 7 b, 43.5 to refer the emission to a 4375 
beam. We divide the SPIRE spectrum by f]b, 43.5 to refer 
the flux density at all wavelengths to that observed by 
a 4375 beam. We also use these values to refer CO in¬ 
tegrated flux values measured from other faci litie s with 
smaller beam sizes to the 4375 beam (Section 2.4). 

We apply an additional correction (also used in K14) 
to match the total flux density of the spectrometer 
with the photometer flux density. At the high fre¬ 
quency end, we match the total SSW flux density inte¬ 
grated over the photometer PSW bandpass, F{PSW), 
to that of the PSW photometer-integrated flux den¬ 
sity at 4375, F"{PSW) by multiplying the spectrum by 

-Assw = F'"(PSW)/F(PSW). There are two photometer 
bands (PMW and PLW) which overlap with the SLW 
band, so we define a line that connects those two ra¬ 
tios, F"(PMW)/E(PMW) and F"(PLW)/ F(PLW), and 
multiply the spectrum by that ratio as a function of wave¬ 
length, AsLw(i^)- This photometry correction step is of¬ 
ten most signiheant in the SSW, which can overestimate 
the measured flux compared to the photometry. Some 
spectra had somewhat over-subtracted telescope emis¬ 
sion, giving slightly negative flux densities, especially at 
the lowest-frequency end. In these cases, no correction 
was done to the SLW band to match the PLW photom¬ 
etry, which would use negative ratios. In a few cases, we 
also did not correct the SLW band if the ratios derived 
from the PLW and PMW were significantly discrepant 
(i.e. would produce a non-sensical SLW continuum). The 
spectra that were not corrected are marked in Table 

In order to compare the CO emission to LfiRj we 
must also correct Lfir to properly represent the same 
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amount of emission as the CO within our beam. Similar 
to the procedure above, we convolve the SPIRE pho¬ 
tometer maps at wavelength A to the beamsize of 43"5 
and find the ratio of the peak flux density in Jy mea¬ 
sured with such a beam (Pbeam.A) to the total integrated 
emission in the map Assuming LpiR^beam = 

ApiR, total ^Ebeam,250/im/Etotai,250/im5 and likewise for the 
350 and 500 maps, we can determine the proper 
ApiR^beam for Comparison to the CO emission. The three 
photometers agreed well, and so we use the average of 
the Ebeam.A/Atotai.A ratios. Of the 232 observations with 
known redshift and available spectra, 118 have ratios of 
< 0.8, and 42 have ratios of < 0.5. We propagate the 
errors from the total measured integrated flux density 
through to the final measurement of TpiR.beam- 

2.3. Herschel FTS Line Fitting Procedure 

The CO J = 4—3 to J = 13—12 lines, both [Cl] lines, 
and the [Nil] 205 /im line are the brightest lines in the 
FTS spectra. To fit these, we start with the FTFitter 
code from the University of Lethbridge]^ Treating each 
detector (SLW and SSW) separately, the code fits a poly¬ 
nomial to the baseline, and then simultaneously fits unre¬ 
solved lines at the expected frequencies of the lines listed 
above (given known redshifts). We place a lower limit of 
the total area of the line profile to be above 0; we do not 
expect any of these lines to be in absorption. We limit 
the line center to within ± 500 kms“^ of that expected 
from the redshift to allow for uncertainty in the velocity 
scale and physical differences in the gas kinematics. In 
wavenumbers, this is about 0.025 - 0.084 cm“^ over the 
band, compared to the FWHM of the line profile of 0.048 
cm“^. 

We manually inspected the resulting fits to determine 
if any lines were clearly resolved. This is most likely to 
be case for the [Nil] line as velocity resolution is the 
highest at the higher frequencies, and it is much brighter 
than CO lines at similarly high frequencies which may 
be undetected. Resolved lines do not show the same 
characteristic ringing of the sine function; the ringing is 
significantly lower, if not imperceptible, smeared out by 
the effective convolution of the emission line profile and 
instrumental profile. We refit the lines that met this cri¬ 
teria as a Gaussian convolved with the instrumental line 
profile. In this case, the lines are barely resolved, thus 
Gaussian profiles are perfectly adequate (no more de¬ 
tailed velocity profiles can be determined from the FTS). 

The fact that SPIRE utilizes a Fourier Transform Spec¬ 
trometer introduces a special problem in the treatment 
of line fitting. The true measured quantity is the in- 
terferogram, or the interference pattern at the focus as a 
function of optical path difference (OPD) as the mirror of 
the interferometer moves linearly. The spectrum itself is 
the Fourier Transform (FT) of this interferogram, which 
leads to two important consequences: 1) the wavelength 
bins are not truly independent, which many fitting rou¬ 
tines assume, and 2) the resulting noise pattern closely 
resembles the FWHM = 0.048 cm“^ sine function line 
profile. The result is that the errors output by a least- 
squares fitter, like the FTFitter and the built-in HIPE 
Spectrum Fitter routines, may not be an accurate rep- 

® https://www.uleth.ca/phy/naylor/index.php?page= 


resentation of the line flux uncertainty. Moreover, it can 
do an excellent job of fitting a “ripple” in the spectrum 
which, to the observer’s eye, may not be particularly dis¬ 
tinguishable from any other ripple nearby, other than 
that we expect the e.g., CO line to correspond to the 
fitted ripple’s wavelength, and no similarly strong lines 
to be adjacent in the spectrum. 

Though the ideal situation would be to the fit the in¬ 
terferogram itself, this is not a user-accessible option for 
SPIRE data considering the many calibration steps that 
occur in processing after the FT. Instead, we created a 
Bayesian analysis method to determine the probability 
distribution function of the true line flux given the ob¬ 
served line flux, P{ftrue\fobs), which is heavily influenced 
by the correlated noise pattern in the spectrum. The 
noise itself is difficult to accurately characterize, varying 
from observation to observation, and across the bandpass 
of a given observation. Therefore, instead of attempting 
to describe the correlated noise for our entire sample, we 
focus on the area around each individual (unresolved) 
line. 

We describe the procedure briefly here, but show 
a more in-depth example with illustrative figures for 
NGC4388 in Appendix This procedure is not used 
for lines that were manumly identified as resolved, which 
are already high signal/noise. For each line, we input 
sine profile lines of varying amplitudes /true over the re¬ 
gion ± 2 cm from the line center (excluding the area 
immediately around any CO, [Cl], or [Nil] lines) and 
then refit the spectrum. We compare tbe measured in¬ 
tegrated fluxes, fobs to the known input values, ftrue- 
The probability distribution function for our CO line is 
the distribution of input fluxes that produced that par¬ 
ticular measured flux value, a slice of the Pifobs\ftrue) 
two-dimensional distribution that we created. 

For high signal-to-noise line detections, the procedure 
replicates a Gaussian distribution of similar median and 
error (a) as the parameters estimated by the FTFitter. 
This is because a very high amplitude line input, added 
anywhere in the spectrum, will return the same inte¬ 
grated flux value as we input (fobs = ftrue)- However, 
a line with smaller amplitude may add constructively 
or destructively to the underlying ripple pattern of the 
spectrum, returning a higher or lower flux than input. A 
local variation from the detector’s average baseline may 
also influence the final fitted value, which may shift the 
median value of the probability distribution function of 
the line flux. 

The procedure makes the most difference for the high- 
J CO lines; 60% of the 3cr detections of CO J = 13—12 
from least-squares fitting were shown to have > 3a 
uncertainty. For all the CO lines in the SSW band 
(>/upper > 9), this number is 44%. For the CO lines 
in the SLW band, it is only 15%. The numbers are the 
lowest for [Nil] (4%) because it is so bright, and for [Cl] 
J = 2—1 and CO J = 7—6 (5%, 8%) because they lie in 
the lowest noise part of the spectrum and are relatively 
bright. 

The results from this line fitting procedure are shown 
in Table The median, — Icr, and -1-lcr values are de¬ 
rived from the points at which the cumulative distribu¬ 
tion functions (CDF) equal 0.5, 0.159, and 0.841, respec¬ 
tively. If — Icr/median is less than 3, a value for a 3cr 
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upper limit is also shown (where the CDF = 0.997). 

TABLE 2 

Line Fluxes and Uncertainty Ranges from SPIRE FTS 


Galaxy 

Line 

Resolyed? 

Median 

Jy km s~^ 

-la 

Jy km s~^ 

H-lcr 

Jy km s~^ 

3cr Upper Limit 
Jy km 

NGC0023 

GIl-O 


2.57e-|-02 

6.81e-|-01 

5.57e-|-02 

7.98e-|-02 

NGC0023 

GI2-1 


3.82e-|-02 

3.09e-|-02 

4.57e-|-02 


NGC0023 

C04-3 


3.12e-|-02 

7.15e-|-01 

6.10e-|-02 

1.21e-|-03 

NGC0023 

C05-4 


7.38e-|-02 

5.76e-|-02 

9.69e-|-02 


NGC0023 

C06-5 


4.23e-|-02 

3.26e-|-02 

4.88e-|-02 


NGC0023 

C07-6 


3.24e-|-02 

2.42e-|-02 

3.97e-|-02 


NGC0023 

C08-7 


2.49e-|-02 

1.18e-|-02 

3.83e-|-02 

5.71e-|-02 

NGC0023 

C09-8 


1.15e-|-02 

2.51e-|-01 

2.94e-|-02 

5.63e-|-02 

NGC0023 

GOlO-9 


2.22e-|-02 

5.53e-|-01 

3.56e-|-02 

5.28e-|-02 

NGC0023 

G012-11 


7.78e-|-01 

1.81e-|-01 

1.92e-|-02 

3.26e-|-02 

NGC0023 

C013-12 


2.29e-|-02 

9.69e-|-01 

3.55e-|-02 

4.76e-|-02 

NGC0023 

Nil 

X 

2.63e-|-03 

2.46e-|-03 

2.80e-|-03 


NGC34 

GIl-0 


4.96e-|-02 

3.36e-|-02 

6.45e-|-02 


NGC34 

CI2-1 


4.29e-|-02 

3.83e-|-02 

4.72e-|-02 


NGC34 

C05-4 


7.87e-|-02 

6.92e-|-02 

8.77e-|-02 



Note. — Table is published in its entirety in the electronic 
edition of The Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 

2.4. Low-J CO Lines from the Literature and Arizona 
Radio Observatory 

The bandpass of the Herschel FTS starts around the 
CO J = 4—3 line, but the majority of the molecular mass 


in galaxies is cool and populates the lower rotational lev¬ 
els. We complement the line fluxes derived from the FTS 
with the CO J = 1—0, J — 2—1, and J = 3—2 lines avail¬ 
able from ground-based observatories. Many of these 
galaxies have already been studied in the literature, es¬ 
pecially in large CO surveys. 


TABLE 3 

CO J= 1—0 TO J= 3—2 Line Fluxes 


Galaxy 

Jup Reported am 

ac. Unit A^ 12^ JA-u 

a 

Ref 



X kms~^ 




NGG0023 

1 

16.9 


4.2 

K 

141 

24 

138.9 

34.7 

1 

NGG0023 

1 

6.0 


1.2 

K 

374 

55 

169.9 

34.0 

2 

NGG0023 

1 

8.9 


2.2 

K 

190 

45 

184.4 

46.1 

1 

NGG0023 

1 

18.0 


1.3 

K 


33 

234.0 

16.9 

3 

NGG0023 

1 

34.0 

0.4 

6.8 

K 


22 

247.4 

49.6 

4 

NGG0023 

2 

18.8 


4.7 

K 

129 

12 

235.6 

58.9 

1 

NGG0023 

2 

7.4 


1.9 

K 

210 

24 

243.2 

60.8 

1 

NGG0023 

3 

15.3 

1.3 

2.3 

K 

257 

22 

1003.0 

171.8 

5 

NGG34 

1 

17.0 


4.2 

K 

149 

24 

115.3 

28.8 

1 

NGG34 

1 

4.5 

0.3 

0.4 

K 

295 

55 

131.8 

16.0 

6 

NGG34 

1 

6.7 


1.7 

K 

274 

45 

138.4 

34.6 

1 

NGG34 

1 

148.5 

13.5 

29.7 

Jy 


45 

147.6 

32.4 

7 

NGG34 

2 

4.3 


1.1 

K 

271 

24 

116.7 

29.2 

1 

NGG34 

2 

56.5 


14.1 

K 

172 

12 

471.9 

118.0 

1 

NGG34 

3 

7.7 

1.8 

1.2 

K 

168 

22 

404.8 

110.5 

5 


Note. — Table is published in its entirety in the electronic 
edition of The Astrophysical Journal. A portion is shown here 
for guidance regarding its form and content. The first eight 
columns refer specifically to the measurements reported 
in the literature. “Reported” is the reported value in the 
units of the “Units” column. <Jm and ctc refer to measurement 
and calibration error, if separately reported, otherwise calibration 
errors are assumed or contain total error. At, is the FWHM of 
the line, if reported, and Q5 is the FWHM of the beam size. The 
next two columns are the values used in our analysis: 
IAv{^fts) is the flux in Jy km/s, referred to the 43 75 beam, 
and a is the total accompanying error. Refe rences, fl) Albrecht] 


let al.l (|2007ll; (21 ISanders et a hllllQQlIl; f3) lElfhagc et a l (|iUUblt7 
(4}|Ciarcia-B uril lo et al.1l|2Ui2|):~T3} bJVLl' (this work): (b)lM aiolmo 
et al.lJTjjjjTl): (7 )|L5a an et al.lJzOOSll: (8) 12M (this work): (^) luerov l 


iMmA (V)l 

rlTOPgi: ff 


’12) |Mira^ et al. 


ni 


et al. J(|l999[l: (jAl 


aolornmi et al. 11 9U7l): (IbllRapadop oulos et aL|Q2Ui2D ; (lUlli Dvans 
et al.NliiUUPil; t2l)lAalto e t al.i riU^jbl); (l^DlLieecn et al .|(|201 0||; (‘2'Z] 
Kame^tzky et al.|(|2Uli|l: (2jTlapinogliQe^STllB^SIn (24)1 Young 


LSayet e't al.l 20061): (11) iMao et al.l 
\ u \ Ir A.,.,,. 11 1 1 I /TmmlrrTTTTm 


lyyut; 113) 


Young et al. 


" Garay et al.l JiUUij ; 

'p95ll: 


[Evans I 


et al 

vm 


] (|2U1 i|l : (23) [^Tnoglio et al.l (|2Ui2[) ; T2TTTYoung| 




B : ('25)lHlatalo A a . |(|20l U: 

ipoulosS^: Seaquis^jlbudlr: ( 


: (2tj)rLazarett et al. 


bandqvi^ j 1 999 


J2003 l ): [ 32) |Yao et a 

|pt a I M9ni A I- LX.'S dSrhirrn pt. a"T 


IS 


(28~)rSanda vist e t al. 
JjothwelT et al, ' 


.'Tggs 


TgggT 


aailkil iwa et al.hl20i2b: (841 i Bosell i 
ii9nun- Vri Xr UwartiTiwinr 


For some galaxies, we also performed single-dish mea¬ 
surements using the Arizona Radio Observatory (ARO). 
Measurements of the CO J — 1—0 line were conducted 
with the 12m dish on Kitt Peak in May of 2015, and 
those of CO J = 2—1 and J = 3—2 were conducted with 
the Submillimeter Telescope (SMT) located on Mt. Gra¬ 
ham from November 2014 to February 2015. At the 12m, 
we used the ALMA Type 3 mm receiver with two 2 MHz 
backends in series, yielding 2.6 kms“^ channel resolu¬ 
tion and about 670 kms“^ bandwidth. At the SMT, 
we used the ALMA Type 1.3 mm sideband separating 
receiver (for CO J= 2—1) and the 0.8 mm double side¬ 
band receiver (for CO J = 3—2) with the 1 MHz hlter- 
banks in 2IF mode. Most observations were conducted 
with beam switching, except for highly extended sources 
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which required position switching. Pointing and focus 
were checked periodically on planets or bright continuum 
sources. 

Beam efficiency measurements were conducted with 
Venus and Jupite^ us ing the procedure described in 
Schlingman et al. (2011). For CO J=l—0, J = 2—1, and 
J = ' 6 —2 we found Tymb of 0.86-0.89, 0.57-0.62, and 0.58- 
0.65, respectively. The beam sizes are approximately 
55", 32", and 22", respectively. The spectra were re¬ 
duced, baseline subtracted, and converted to Tmb scale 
using rjjnh in CLASS. Finalized spectra were smoothed 
to approximately 20 kms“^ bins from which total inte¬ 
grated fluxes were derived. 

All low-J lines utilized in this work, including the ones 
from ARO, are available in Table As Table shows, 
the low-J lines we use come from a variety of telescopes 
with different beam sizes. For subsequent comparison 
to Herschel CO lines, all line fluxes are referenced to the 
43'.'5 beam size using the same ratios {rjb, 43 . 5 ) as described 
in Section 2.2 Table lists both the literature reported 
values (third column, in their original beam sizes and 
units) and the 4375 referenced fluxes in Jy kms“^ we 
use in our analysis (ninth column). 


3. ANALYSIS 


We examined the relation ship between L^rluand Lftr 
( from 40-1 20 qm), similar to Lu et al. ( 2014|), iGreve et al. 
(20141) andlLiu etal. I (120151), discussed turthi 


er in Section 

3.2 ~We choose the orientation of our axes, with Lfir on 
the y-axis, for the easiest comparison to the existing lit¬ 
erature. This orientation comes from the comparison to 
the Kennicutt-Schmidt (K-S) scaling law, which relates 
star formation rate (SFR) to molecular ga s surface den¬ 
sity, with S* oc (Kennicutt||l998). Subsequent 

large-scale analyses have found indices from 1.0-1.4 for 
the molecular gas and higher values, 1.4-3.1, for the to¬ 
tal gas surface density ( Kennicutt fc Evans||2012 1. Lfir 
can be considered a proxy tor SFR it one excludes the 
contribution of AGN to Lfir, which we admittedly do 
not separate here. It is likely small for most sub-ULIRG 
galaxies, especially at these wavelengths. In the case of 
J = 1—0, the well-known “X-factor” or aco is used to 
relate the luminosity of CO J = 1—0 to the total molecu¬ 
lar mass, so the slope derived here is comparable to the 
K-S relation. However, this relationship is not applicable 
at higher-J, where the CO luminosity is not a tracer of 
mass (K14), but we choose to keep the same orientation 
to avoid confusion. Neither variable should be consid¬ 
ered the “independent” one, regardless of which appears 
on the x-axis. 

The theoretical explanation for these relationships 
were Hrst investigated to describe the discrepant power 
laws in low-J emission of various molecules, namely CO 
J = 1—0, where Lft roc and HCN J = 1—0, 

where Lfiroc Lhcn ( Gao fc Solomon|[2004a|b ). The CO 


^ Additional planetary observations provided by Karin Sand- 
Strom. 

® Lqq= 3.25 xl0'^ScoAnD|(l + [K kms^^pc^], 

where is the In minositv distance in Mn e in GHz. Snr,A v 

in Jy kms“^, from 


ninositv distance in Mn e in GHz. SnnA v 

Carilli fc Walter] | |2013 i; [Solomon et al.|jl992| . 


Note that Lqq is just the area within the beam times the veiocity- 
integrated antenna temperature, A X TAv, where A is the area 
10l2r2D2(l + 2)-4 [pc2]. 


slope closely resembles that of the aforeme ntioned K- 

5 relation slope, but HCN does not match. Krumholz 

6 Thompson (2007) showed that this could be under¬ 
stood as a consequence of the different critical densities 
for different species’ ground-state transitions, assuming 
isothermality. In short, a molecular line with a low crit¬ 
ical density compared to a galaxy’s median gas density, 
Pg, will be excited by the majority of the molecular gas. 
The star formation rate will therefore depend on one fac¬ 
tor of density pg based on the total amount available 
for star formation, and a factor from the dynami¬ 
cal timescale of the gas, adding up to a total factor of 
1.5 in the case of the low-ncrit CO J = 1-0 line. On the 
other hand, when the molecular line has a critical density 
higher than that of the median gas density, its emission 
will be picked out from high-density peaks only, specih- 
cally peaks of the same density (and therefore same free- 
fall time) across different galaxies. Therefore the higher 
ricrit of HCN J = 1-0 yields a slope of 1.0. We discuss 
our results in context of these expectations in Section 

3.1. Fitting Lqq/Lfir Slopes 

As mentioned in Section [272} all fluxes including Lfir 
are referred to the emission within a 4375 beam. To 
determine the coefficients of the relation log(LFiR) = a 
log(LQQ) + b, we used the python module Inr.bces]^ 
which utilizes the Bivariat e, Correlated Errors and in - 
trinstic Scatter method of Akritas & Bershady ( 1996|). 
This is important because we have errors on both vari- 
ables (we introduced a non-negligible error into the Lfir 
variable through our source/beam correction). As stated 
above, we chose the examination of Lfir as a function of 
^CO to match the most recent literature and note that 
the solution to the inverse problem (Lqq as a function 
of Lfir) does not simply produce best-fit slopes that are 
the inverse of those presented here. For the case of low-J 
lines collected from the literature and our ARO follow¬ 
up, multiple measurements for the same galaxies (posi¬ 
tioned at the location matching the FTS coordinates) are 
treated independently. 


TABLE 4 

Correlations between and Lfir: Full Sample. 



a 

CTa 

b 


n 

1 

1.27 

0.04 

-1.0 

0.4 

299 

2 

1.11 

0.07 

0.6 

0.7 

138 

3 

1.18 

0.03 

0.1 

0.3 

131 

4 

1.09 

0.05 

1.2 

0.4 

108 

5 

1.05 

0.03 

1.8 

0.3 

195 

6 

1.04 

0.03 

2.2 

0.2 

199 

7 

0.98 

0.03 

2.9 

0.2 

196 

8 

1.00 

0.03 

3.0 

0.3 

186 

9 

1.03 

0.04 

2.9 

0.3 

176 

10 

1.01 

0.03 

3.2 

0.3 

184 

11 

1.06 

0.04 

3.1 

0.3 

166 

12 

0.99 

0.03 

3.7 

0.2 

168 

13 

1.12 

0.04 

2.9 

0.3 

156 


Note. — Best fit measurements and errors for Log(LFipi)=a 
Log(Z;QQ)-|-b. Column n = number of data points used in relation 
(not all 3cr detections). 


® http://home.strw.leidenuniv.nl/-sifon/pycorner/bces/ 
by Cristobal Siton 
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log L' [K km/s pc" ] log L' [K km/s pc^ ] log L' [K km/s pc" ] 


Fig. 1.— CO vs. Lfir- The y-axis is the Tpir in the beam for comparison to the CO measurement (see Section |2.2| Low-J lines 
may include multiple measurements for the same galaxy if available in the literature (Table [^. Blue data points indicate resolved (Sinc- 
Gaussian) measurements from the FTS. Line fits are described in Section]^ The green line is fitting the whole sample; the CO J = 1—0 
line fit is shown as a dotted black line on each other CO plot for comparison. The dashed orange and purple lines are fits when separating 
the sample into galaxies above and below Lfir= 6 X 10^*^ Lq , respectively. 


When including all spectra in our sample, we find 
slopes starting at 1.3 for CO J = 1—0, and lowering to 
about ~ 1 for the mid- to high-J CO lines, with no dis¬ 
cernible trend with increasing J. The results are shown 
in Table There is no significant difference whether we 
include or exclude FTS lines with S/N < 3. We also sep¬ 
arated our samples into a few categories that are some¬ 
what overlapping, and summarize the differences here: 


We separated our galaxies into known AGN (cat¬ 
egorized in Hyperledap^ 
of Seyfert) or not. This 


as quasar or any type 
classihcation (using the 


agnclass catego ry) is taken from the |V eron-Cetty 


& Veron (|2006|) catalog. Within those classihed 
by Hyperleda as AGN (78 of the 232 galaxies), no 
information is provided regarding the relative SF 
vs. AGN contributions to the total Lfir, so this 
is a somewhat crude division of the galaxy sam¬ 
ple. Looking only at AGN (Table [^, compared to 
the whole sample we find higher slopes for the low- 
J lines (1.5 ± 0.1, 1.2 ± 0.1, 1.3 ± 0.1, 1.1 ± 0.1, 
1.15 ± 0.04 for J = 1—0 through J = 5—4), but the 
slope error bars overlap by CO J = 6—5 and con¬ 
tinue to follow the same trend as the whole sam¬ 
ple. Looking at the sample that completely ex¬ 
cludes AGN, we find a lower slope than the whole 


http://leda.univ-lyonl.fr/leda/rawcat/al09.html/ 


sample for CO J=l—0 (1.13 ± 0.06), but at subse¬ 
quently higher-J the slopes are not distinguishable 
from the combined sample. In summary, compar¬ 
ing AGN to the non-AGN sample, the most signifi¬ 
cant difference is in the GO J = l—0 line. There are 
also differences, at less significance, up to the J = 
5—4 line. The low-J CO emission is not expected 
to be affected by the AGN; it is the Lfir and high- 
J CO likely being influenced. We likely do not see 
any difference because our AGN-designated galax¬ 
ies (which are only about one third of the full sam¬ 
ple) may not be entirely dominated in their molec¬ 
ular excitation from the AGN; as mentioned, we 
do not separate the relative SF vs. AGN contri¬ 
butions to total Lfir. Better quantification of the 
AGN influence, and higher spatial resolution, may 
result in differences in the slope. 


Astronomers often separate galaxies into (U)LIRGs 
or lower luminosity galaxies. We separate our 
galaxies at Lfir = 6 x 10^° Lq which we found is 
approximately equal to Lir (8-1000 /r m) = lx 10^^, 
based on the luminosities listed in IGreve et al.l 


(2014) and Kamenetzky et al. 
cutoff value does not change 
follow.) The CO 


(2014). (The exact 


the conclusions that 
J = 1—0 line has been known to 
not be fit by one slope among (U)LIRGs and lower 
luminosity galaxies; the superlinear slope that re¬ 
sults from a single power line ht is due to higher 
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dense gas fractions in (U)LIRGs 1 

Gao & Solomon 

2004b 

Greve et al. 

2014|), esseni 

tially creating a 


slopes in J= 1—0 among these two populations fit 
separately than their combined slope of 1.3 ± 0.4. 
We find the mid-J lines of CO are well fit by a single 
slope across many orders of magnitude; the slopes 
in all three cases (full sample, just (U)LIRGs, just 
galaxies with lower luminosities than (U)LIRGs) 
are not statistically distinguishable given the error 
bars. Our three highest-J CO lines (J= 11—10, J = 
12—11 J = 13—12), however, do show a measurable 
difference in the best fit slope. Focusing primarily 
on (U)LIRGs changes the slopes at higher-J, de¬ 
creasing to about 0.83 ± 0.03 (weighted average of 
three highest-J CO lines in Table , a highly sig¬ 
nificant difference from 1. The results for these two 
populations are shown in Tables and 


• We restricted the fit to only those galaxies not par¬ 
ticularly well resolved, those with TFiR,beam > 0.8 
J^FiR,total based on the photometry maps, to see if 
our infrared luminosity correction could be influ¬ 
encing the y-axis values. We find lower slopes in 
this case, 0.88 ± 0.06 on average for high-J lines, 
but this may be due to the fact that 2/3 of this 
population are (U)LIRGs, not because large Tfir 
corrections are necessarily inaccurate. 


Certainly, these slopes are masking a considerable 
amount of scatter in the data, and different trends can 
b e discerned when fitting different populations (also seen 
Liu et al. (2015|)). The remaining figures of this paper 


examine the survey populations themselves with compar¬ 
isons to these derived, “global” relations. 


TABLE 5 

Correlations between I,/,q and Lpir: AGN Only. 


J up 

a 

CJ-a 

b 

O-b 

n 

1 

1.54 

0.07 

-3.4 

0.7 

99 

2 

1.24 

0.14 

-0.6 

1.2 

39 

3 

1.30 

0.08 

-1.1 

0.7 

39 

4 

1.15 

0.08 

0.6 

0.6 

33 

5 

1.15 

0.04 

0.9 

0.3 

60 

6 

1.09 

0.04 

1.6 

0.4 

62 

7 

1.03 

0.04 

2.5 

0.3 

61 

8 

1.03 

0.05 

2.6 

0.4 

58 

9 

1.02 

0.06 

2.9 

0.5 

53 

10 

1.01 

0.05 

3.1 

0.4 

58 

11 

1.01 

0.07 

3.3 

0.5 

56 

12 

0.99 

0.06 

3.6 

0.4 

55 

13 

1.13 

0.07 

2.7 

0.5 

53 


Note. — Best fit measurements and errors for Log(LFiR)=a 
Log(LQQ)-|-b. Column n = number of data points used in relation 
(not all 3 (t detections). 


TABLE 6 

Correlations between AND Lpir: (U)LIRGs Only. 


«^up 

a 

CTa 

b 

O-b 

n 

1 

1.15 

0.09 

0.2 

0.8 

225 

2 

0.66 

0.15 

4.9 

1.4 

100 

3 

0.94 

0.11 

2.3 

1.1 

86 

4 

0.68 

0.13 

4.9 

1.1 

44 

5 

0.96 

0.07 

2.7 

0.6 

129 


TABLE 6 — Continued 



a 

O'a 

'b 

<^b 

n 

6 

1.02 

0.05 

2.4 

0.4 

132 

7 

0.92 

0.04 

3.5 

0.4 

131 

8 

0.91 

0.04 

3.7 

0.3 

125 

9 

0.92 

0.07 

3.8 

0.5 

111 

10 

0.83 

0.03 

4.7 

0.3 

126 

11 

0.86 

0.05 

4.7 

0.4 

115 

12 

0.79 

0.04 

5.3 

0.3 

114 

13 

0.85 

0.05 

5.0 

0.4 

109 


Note. — Best fit measurements and errors for Log(LpiR)=a 
Log(Z/QQ)-|-b. Column n = number of data points used in relation 
(not all 3cr detections). 


TABLE 7 

Correlations between L/,q and Lpir: Non-(U)LIRGs Only. 


'^up 

a 

O-a 

b 

^b 

n 

1 

1.05 

0.09 

0.8 

0.8 

74 

2 

1.12 

0.11 

0.4 

0.9 

38 

3 

1.05 

0.05 

0.9 

0.4 

45 

4 

1.09 

0.07 

1.2 

0.5 

64 

5 

1.01 

0.06 

2.1 

0.5 

66 

6 

1.01 

0.05 

2.3 

0.4 

67 

7 

1.00 

0.06 

2.8 

0.4 

65 

8 

1.07 

0.07 

2.5 

0.5 

61 

9 

1.06 

0.08 

2.7 

0.5 

65 

10 

1.12 

0.07 

2.5 

0.5 

58 

11 

1.10 

0.09 

2.7 

0.6 

51 

12 

1.03 

0.07 

3.4 

0.4 

54 

13 

1.23 

0.08 

2.2 

0.5 

47 


Note. — Best fit measurements and errors for Log(LpiR)=a 
Log(LQQ)-|-b. Column n = number of data points used in relation 
(not all 3cr detections). 

Fitting a linear relation in log space requires creating 
error bars that are symmetric in log space, which ours are 
not. The exaggeration of the linear errors when viewed 
in log space is highest for those line measurements with 
the lowest signal-to-noise ratio. To test whether the con¬ 
version to these symmetric error bars (which decreases 
the size of the lower error bar) introduces systematic bias 
in deriving the slopes, we chose fixed values of the slopes, 
set the nominal y-values to match the chosen slope, and 
draw randomly from each data point’s x and y error bars 
in linear space. Fitting such scattered distributions many 
times results in distributions centered upon the chosen 
slopes and with scatter comparable to our error bars in 
Table We find no evidence of systematic bias due to 
this effect. 

We also investigated whether our error bars are respon¬ 
sible for the slightly lower value of the CO J = 1-0 line’s 
slope at 1.3 instead of 1.4-1.6. This value is not due to 
the error bars in either the x- or y-directions; we still 
find a slope of about 1.3 if we exclude one, the other, or 
both in the line fitting. We do find a slope of 1.44 if we 
restrict the fit to those data points with log {L'qq) > 9, 
indicating that the lower luminosity points are bringing 
down the slope overall. 


3.2. Comparison to Similar Works 


Greve et al. (2014) conducted a comparison of CO 
emission to Lftr tor the non-extended galaxies o f the 
HerCULES sample (|van der Werf et al.| j2010|) and 
found sharply decreasing slopes between log(LFiR) and 






































Lqq/Lfir Relations Across the Herschel SPIRE Archive 


11 


log(LQQ), starting around 1 at CO J=l—0 and decreas¬ 
ing as one moves to higher-J starting around CO J = 
6-5 (0.93 ± 0.05, down to 0.47 ±0.20 by the J= 13-12 
line). 

They interpret this sublinear slope as an indication 
that FUV radiation fields are not responsible for the the 
dense, warm gas emitting in the high-J lines. Mechani¬ 
cal heating and shocks are the likely explanation for the 
high-J excitation, as was shown in many galaxies of the 
sample of K14 (and references therein). We are in agree¬ 
ment with Greve et al. (2014) that mechanical feedback 
can be related to star formation, such as stellar outflows, 
winds, and supernova remnant expansion; not being re¬ 
lated to the FUV radiation from stars does not mean the 
excitation is not relat ed to star formation at all. 

The fluxes used in | Grev e et al. (2014) are reported 

I'heir sample is included 


J CO slopes of 1 in a large sample of galaxies and less 
t han 1 for (U)LI RGs. 


Rosenberg et al. (2015). 


m 

in ours, but our measured fluxes do not match in all 
cases. For lines up to J = 10—9, 2/3 of their 18 non- 
extended galaxies match our flux values within 30%. At 
higher-J, this number is more like 1/3. When discrepant, 
their values are often ^ 50% higher, or even factors of a 
few for IRASF05189-2524, for which we do not use the 
same spectra. The major difference for other galaxies is 
likely the treatment of t he sou rce/beam correction. Im¬ 
portantly, |Greve et al.| (|2014|) only included (U)LIRGs 
(log(LFiR) > 11)- As discussed above, we find lower 
slopes if we restrict ourselves to this population, bu t not 
as low as the slopes reported by Greve et al. (2014), due 
to the differences in our measured fluxes. 


Liu et al. (2015) conducted a larger study, more com- 
parable to ours, utilizing the full extent of the Herschel 
archive. With the inclusion of galaxies from log(LFiR) 
> 8, they found linear slopes between log(LFiR) and 
log(LQo) throughout the CO ladder starting at J = 4—3 
(statistically consiste nt with our results) . Also consistent 
with our results and Greve et al.| (2014), restricting the 
fits to the HerCULES sample of (U)LiRGs yields sublin¬ 
ear slopes for the highest-J lines. Their work does not 
include a table of the galaxies included in their sample 
or the line fluxes used, so we cannot make a direct com¬ 
parison of the fluxes used for our relations at this time. 
Some major differences may be important in the com¬ 
parison of our work: first, they use SPIRE calibration 
version 12 and we use version 13, which may make the 
most difference for galaxies observed near the start of 
SPIRE cooling cycles (due to “cooler burp”). Second, 
each of their GO and Tfir relations use a different beam 
size, that of the frequency of the CO line in the FTS. 
Third, for sparsely and intermediately sampled galax¬ 
ies, they extract spectra from multiple bolometers in the 
detector array, meaning they are including multiple spa¬ 
tially Nyquist-sampled data points in their relationships 
(with their own matching Lfir values) for such galaxies. 
Finally, their line fitting procedure uses sinc-convolved- 
gaussian (SCG) lines in HIPE for sparse observations and 
sine for intermediate or full sampling observations. They 
note that SCG-derived fluxes are systematically higher, 
but we found that few galaxies have CO linewidths large 
enough to be resolved by the FTS. Diffe r ences in the 


slopes of (U)LIRGs among Liu et al. (2015), Greve et al. 
(2014), and this work can also be attributed to the small 


dynamical range spanned by (U)LIRGs. Despite these 
differences, our results agree well in finding mid- to high- 


Lu et al. (2014) also examined the GO J = 4—3 through 
J= 13—12 emission of the 65 LIRGs in the Great Obser¬ 
vatories All-Sky LIRG Survey (GOALS), comparing to 
the IRAS 60-to-100 ^m color, C(60/100), as a proxy for 
dust temperature. They demonstrated that Lir is not 
the best predictor of SLED shapes; we find an overall 
trend in Figure (the line luminosity ratios relative to 
J = 1-0 in average SLEDs by Lfir bin increase with 
^firJi but also a significant amount of variation in Fig¬ 
ure ^ (individual mid- to high-J CO/J= 1-0 luminosity 
ratios for each galaxy, unbinned), discussed in the next 
section. 0(60/100), which is not presented here, is po¬ 
tentially a better predictor of the GO SLED shape (based 
on the line at which the SLED peaks in luminosity). 

3.3. Average SLEDs by Lfir Ranges 

Within galaxy-wide log(LFiR) ranges of approx¬ 
imately 0.5 dex, we compiled weighted-averaged 
{Leo/L fir)/{L co{.J = 1-0)/Lfir) ratios, presented in 
Figure (top; see caption for note about conversion to 
SLEDs are normalized to the J = 1—O line to 
show the relative excitation across the CO ladder. The 
CO J = 1-0 line measures emission from the same type 
of cold, ubiquitous molecular gas found throughout many 
types of galaxies. There are four main trends to examine 
in this plot: First, there is a trend with increasing Lfir 
towards much more high-J CO luminosity compared to 
J = 1—0. This indicates a greater energy input relative 
to typical PDRs, to explain the high-J emission. Second, 
the location of the CO luminosity peak moves to higher- 
J with higher Lfir- Third, the slope of the mid- to 
high-J CO emission SLED becomes flatter with increas¬ 
ing Lfir- However, all the SLED slopes are relatively 
flat, none show an extreme drop-off. Fourth, the values 
of mid- to high-J CO relative to CO J = I—0 in lumi¬ 
nosity only range from about 10 to 100 across all lines p] 
This is consistent with the CO SLED compilation shown 
in Figure 5 of K14, which showed high-J/J = 1-0 ra¬ 
tios from about 5 to 100 (with one outlier, NGC6240, at 
250 at its peak, see Meijerink et al. (2013)). These last 
three trends are indicative of a higher average pressure 
(product of kinetic temperature and density) required to 
explain the shape of the CO SLED. We emphasize av¬ 
erage because all of the CO emission is the sum of a 
gradient of conditions from different environments. 

Figure (top) also shows a difference between the 
bin-averaged SLEDs (with data points) and the SLEDs 
that would be predicted simply from our log(LFiR) vs. 
log(L(;;o) slopes (light colored, no data points). The 
predicted SLEDs span a narrower range than the bin- 
averaged SLEDs. Describing each CO line by a single re¬ 
lationship over the entire sample averages over the very 
real diff eren ces in populations, such as those described in 
Section 3.1 (e.g. different Lco/Lfik slopes at higher lu¬ 
minosities]. Also, the actual weighted averages are often 
influenced by high S/N galaxies that may not represent 
the whole bin, but still illustrate the overall trend. Ei¬ 
ther way of examining the data, which we do more in 

In brightness temperature units (L(-,q), this is equivalent 
to 0.16-1.6 for J = 4—3 and 0.0046-0.046 for J = 13—12, because 

Lqq. 
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Fig. 2.— Top: average SLEDs by Lfir ranges. For Lfir bin 
ranges shown in the legend, the value of Lqq /Lyir was averaged if 
measurements existed for at least 3 galaxies. All SLEDs were then 
divided by the value of the Lq,o(J = / Lrir line to demonstrate 

the difference in relative excitation (shape) of the CO ladder. The 
number of data points used in each SLED may change with each 
line, which is why a range is given in the legend. The highest lu¬ 
minosity bin is dominated by more distant galaxies, where the CO 
J = 4—3 line is likely to be redshifted out of the FTS band, which 
is why that black data point is missing. The lighter lines with no 
data markers indicate the predictions from the slopes in Table[^for 
the center of each log bin. Bottom: comp arison to theoretical mod¬ 
els. T he red model predictions are from [Narayanan Krumhol^ 
( |2014|| , for SFR surface densities of 1 (dashed) and 100 (solid) M© 
kpc~^ vr~^ (up to J=_9— 8). The black model predictions are from 
[Kazandjian et ah] l |2015[ | for solar metallicity and Ay = 10. The 
first two predictions are for ngas = lO'^, Go = 10^, o = 0 (dash- 
dot-dot, which drops quickly) and a = 0.1 (dashed, which rises too 
high). The solid black line the sum of two models, rigas = 10^, 
Go = 10^, Q; = 0 to fit the lowest-J lines, and then the model of 
rigas = 103, Go = 10®, a = 0.5, multiplied by 0.1, to attempt to 
(though not well) reproduce a flat mid- to high-J spectrum at less 
than 100. The lines with data points correspond to the average 
SLEDs in the top panel. Note: Lqqoc Lqqoc vIco [Jy 

kms“^]. 

the following section, shows a relatively narrow range of 
high-J to CO J = 1—0 ratios compared to expectations 
from some theoretical models. 

4. DISCUSSION 

In the previous section, we explained the observa¬ 
tional and theoretical motivations for fitting the slopes of 
log(LFiR) vs. log(Lco)- We now place the trends uncov¬ 


ered in the previous section in the context of theoretical 
predictions for CO emission. One may initially attempt 
to explain our linear slopes for the mid- to high-J CO 
lines using the same logic as already presented in Sec¬ 
tion^ to describe the linear slope of the HCN J = 1—0 
line (because the critical densities increase with upper-J 
level). This cannot naively be touted as the explanation 
for our linear slopes because the assumption of iso ther- 
mality is not correct (Krumholz & Thompson|2007 Sec¬ 
tion 4.3.2). Their models rely on the gas temperature 
being lower than the upper-state energy level of the tran¬ 
sition. The modeling of K14 has already shown that the 
kinetic temperature of the high-J emitting gas is quite 
high and that the lines are not thermalized; both temper¬ 
ature and density play an important role in the emission. 
The conditions are not uniform, and this sensitivity to 
both temperature and density (and invalidity of LTE) 
is why the high-J CO emission is not linearly related 
to warm molecular mass. Such a relationship (between 
high-J CO emission and mass) was not found in K14 for 
those reasons. Excitation modeling assuming statistical 
equilibrium parameterized by pressure (temperature x 
density) can better illuminate the physical conditions of 
t he molecular gas. 

[Narayanan et ah (|2008|) ext ended the argument of 
Krumholz & 'Thompson ”(|2007 ) by applying excitation 


and radiative transfer calculations to hydrodynamical 
simulations of disk and merger galaxies. They matched 
the known relations and predicted those for the (at the 
time) unobserved mid-J CO to be less than linear, de¬ 
creasing from about 0.6 to 0.4 from J = 4—3 to J=7—6, 
based on the density and temperature distributions de¬ 
rived from their simulations. For these lines, the decrease 
in the slope is due to the fraction of the emission dom¬ 
inated by subthermal excitation increasing with critical 
density of the tracer. This reiterates the above point 
that any one individual high-J CO line is a poor tracer 
of mass. We do not find such low slopes, and also do 
not discern any trend after J = 5—4. In these relations, 
they strictly considered the SFR, not Tfir (which we 
a nd others use as a proxy). 


[Narayanan & Krumholz (2014) showed that the CO 
SLED shape in their models can be parameterized by 
star formation rate surface density. Their simulations 
only considered heating by far ultraviolet (FUV) pho¬ 
tons, cosmic rays, and energy exchange with dust. FUV 
heating is the driving force behind photodissociation re¬ 
gions (PDRs). Two examples are shown in Figure]^ (bot¬ 
tom); some star formation rate surface densities are a 
qualitative match to our highest luminosity (and there¬ 
fore highest SFR) galaxies, but only up through the mid- 
J lines. Though they do not predict above J=8—7, the 
models imply a sharp drop-off would begin above this 
line. However, even in galaxies of two orders of magni¬ 
tude lower luminosity, we still find quite flat SLED slopes 
in the high-J lines that are not well matched by the mod¬ 
els, which drop off steeply after mid-J. 

As already mentioned in the introduction, typical 
PDRs cannot account for the highly excited CO SLEDs 
seen in many Herschel spectra. The need for energy 
sources beyond FUV photons wa s demonstrated with 
data available from the ground (jPapadopoulos et al.j 


2012), but Herschel has made it even mor e clear in 
a number of individual studies and surveys ( Kamenet-| 
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10 ® 10 ® 10 ^° 10 ^^ 10 ^^ 10 ® 10 ® 10 ^° 10 ^^ 10 ^^ 10 ® 10 ® 10 ^° 10 ^^ 10 ^^ 10 ® 10 ® 10 ^° 10 ^^ 10 ^^ 


LpiR [L„] in Beam 


Lpij, [L„] in Beam 


Lpip [L„] in Beam 


Lpij, [L„] in Beam 


Fig. 3.— -^'Co/^CO,!—0 -^fir- The x-axis is the LpiR in the beam for comparison to the CO measurement (see Sectionj2.2| The 
y-axis is the luminosity ratio compared to J= 1—0 for each line. Black data points indicate 3cr detections in both lines; gray indicate less 
than 3 sigma in the higher-J line. The dotted line denotes Jupper^ the theoretical level for thermalized emission (off the top of the y-axis 
after J = 6—5). The solid line denotes the ratios based on the line fits in Table[^(not a fit to these data, and not utilizing the same data; 
this plot relies on 1:1 matching of CO lines with a J= 1—0 measurement from the same galaxy, and multiple low-J lines from the literature 
are averaged together). Line ratios across the range of Lfir remain generally the same across lines above J=6—5, consistently between 1 
or 10 to 100. 

minosity relative to the J = 1—0 line. Two examples (out 
of a much larger parameter space) are shown in Figure 


zky et al. 


|2014[ Lu et al!]|2014[ [Pereira-Santaella et al. 


201^ and individual references within all three). IT 


et al. (|2014|) divided the shock scenario into two cafe^ 
gories: those associated with current star formation (su¬ 
pernovae, winds) and those associated with other non- 
SF-related phenomena (AGN-driven outflows, radio jets, 
or galaxy-galaxy collisions). They found SF-dominated 
galaxies all had similar ratios of total mid-J {J = 5—4, 
J=6—5, J=7—6, J = 8—7, and J=10—9) CO luminos¬ 
ity to Air, whereas galaxies with non-SF-related shocks 
and high AGN contribution had higher and lower ratios, 
r espectively. 


I 


r 


iKazandjian et al. (20151 extended the treatment of 


PDR models to include varying degrees of influence 
from mechanical heating and predicted the CO emission. 
They found that high-J CO line ratios are especially sen¬ 
sitive to mechanical heating; the same cloud conditions 
(parameterized as the gas density and FUV irradiating 
flux, n and Gq) can produce dramatically different CO 
SLEDs with only a few percent of the total heating con¬ 
tributed by mechanical energy (a = Tmechanicai/Ptotai) 
such as turbulence and winds (see their Figure 6). We 
examined the CO high-J/ J=l—0 luminosity ratio s and 
attempted to compare the|Kazandjian et al. (20151 mod¬ 
els to our overall trends. Vvhile the addition of mechani- 
cal energy {a > 0) can result in flatter high-J spectra, we 
find these ratios dramatically overpredict the high-J lu- 


(bottom) to illustrate the impact of a. With no me- 

lanical heating, the shape of the SLED falls down too 
dramatically at mid-J. With the addition of mechanical 
heating, though the shape is flatter, it rises far above the 
ratios we find (and off the top of the plot), which gener¬ 
ally do not go above 100 Aco,i-o, and never above 180 
Leo, 1-0 for 3(7 detections. We can somewhat reproduce 
the average shape by combining multiple models with 
the higher excitation component reduced by a large per¬ 
centage, indicating a negligible amount of CO J = 1-0 
emission from this component. One example is shown in 
Figure but more detailed comparisons to the models 
examining all the possible parameter space of the models 
is required. 

Examination of Figure [^indicates that our trends de¬ 
rived from the Lrir/Lqq fitting (solid, light colored 
lines) may not be representative of the population as a 
whole, given the variance in the average SLEDs from 
these trends. We therefore also examined the line ratios 
Lco/Lco,i-o as a function of Lrir in FDurej^ These 
data are different than that fit in Figure [^because a) it 
relies on the high-J lines for a given spectrum having a 
corresponding CO J = 1-0 line and b) it averages the 
CO J = 1-0 luminosity when multiple measurements are 
available (though all are referenced to the same beam 
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size). The same behavior we see in the average SLEDs is 
still present; the luminosity of the high-J lines relative to 
CO J = 1—0 only varies by about 1.7 orders of magnitude 
across the range of Lfir- This is true for each of the CO 
lines, which is why the slopes in the average SLEDs are 
all relatively flat. This way of looking at the data also 
illustrates the differences in population, e.g. (U)LIRGs 
lying above the average trendline, which are averaged out 
when fitting the log(LFiR)/log(Lco) slopes as in Figure 

m 


We used the same method as in K14 to conduct 2- 
component non-LTE likelihood modeling of the average 
CO SLEDs by the same Lfir bins as in Figure!^ We de¬ 
scribe the SLED as a sum of two components ofgas, each 
described with 4 parameters: kinetic temperature, den¬ 
sity of molecular hydrogen, column density of CO, and 
area filling factor. While the molecular gas exists in a 
continuum of conditions, two components is the simplest 
description of these conditions (see further discussion in 
K14). The product of temperature and density, the pres¬ 
sure P/k in K cm“^, largely determines the relative shape 
of the SLED. The product of column density and area 
is proportional to the mass, which determines the total 
emission (and as previously discussed, in the case of the 
cold gas only, is often considered directly proportional to 
CO J=l-0). 

The trends with Lfir for the predictions from the 
Lco/Lfir slopes (light solid colors in Figure can 
be mostly foreseen from the shapes alone; the pressure, 
mass, and luminosity for both the warm and cold compo¬ 
nents increase slightly with increasing Lfir, but overall 
the physical conditions are not too dissimilar up through 
log(LFiR) = 12- The ratio of warm to cold CO lumi¬ 
nosity is about 45-65 for all bins log(LFiR) < 12, but 
higher for ULIRGs, around 200. The warm/cold compo¬ 
nent pressure on average drops from about 60 to 50 from 
the log(LFiR) = 9.5-10 bin to the log(LFiR) = 11.5-12.0 
bin, but is only about 25 for the highest bin of log(LFiR) 
> 12. The cold component pressures range from log(P/k 
[K cm—3]) = 4 to 4.5, but 4.7 for ULIRGs, and the warm 
component from 5.7 - 6. Within this parameter, there is 
substantial degeneracy between temperature and density. 
Finally, the warm/cold mass ratio is about 0.2, but only 
0.1 for the ULIRG bin. This means that while both the 
cold and warm components have higher pressure overall, 
more of the mid-J emission in the ULIRGs can be ex¬ 
plained by higher bulk excitation of the total molecular 
gas (most of which is cold). This could have implica¬ 
tions for the stellar initial mass function (IMF) and its 
subsequent effect on the surrounding gas. 

In summary, we find linear slopes between mid- to 
high-J log(LgQ) and log(LFiR). Because this CO emis¬ 
sion is not thermalized, it should not be used as a proxy 
for molecular mass or a strict K-S relation. Such slopes 
are inconsistent with th e hydrodynamical simulations of 


Narayanan et al. (2008). The relative luminosity of high- 
J CO to J = 1-0 slightly correlates with LfiR) but only 
ranges from about 10 to 100 for J = 6—5 through J = 
13—12 when a single power law describes each line. When 
examining the full sample (Figure]^, this range varies 
from a few to 180. Across the range of Lfir here, the 
SLEDs relative to CO J = 1-0 are fairly flat across these 
lines. Neither hydrodynamical nor PDR-|-mechanical 


heating models reproduce this trend when used as a sin¬ 
gle descriptor of the galaxy-wide emission. Combinations 
of such descriptors, which essentially adjusts the relative 
contributions of molecular gas components, could better 
describe the SLEDs. 

5. CONCLUSIONS 

We have presented a catalog of all CO, [Cl], and [Nil] 
lines available from extragalactic spectra from the Her- 
schel SPIRE Fourier Transform Spectrometer. 

1. Our catalog includes a uniform treatment for 
source/beam coupling correction and proper es¬ 
timates of the probability distribution functions 
(PDFs) of line flux measurements given the highly 
correlated nature of the Fourier Transform Spec¬ 
trometer. 

2. Relations of the form log(LFiR) = a log(LQQ) -|- 
b show linear slopes over multiple orders of mag¬ 
nitude for mid- to high-J CO lines, and slightly 
sublinear slopes if restricted to (U)LIRGs. 

3. Average SLEDs show increasing mid- to high-J CO 
luminosity relative to CO J = 1—0, from a few to 
^ 100, with increasing Lfir- Even for the most lu¬ 
minous local galaxies, the high-J to J= 1-0 ratios 
do not exceed 180. 

4. The luminosities relative to CO J = 1-0 remain 
relatively flat from J = 6—5 through J = 13—12, 
across many orders of magnitude of Lfir- 

5. Current theoretical models do not reproduce such 
flat SLEDs with ratios < 180 across such a large 
range of galaxy luminosity. 

6. Preliminary RADEX modeling shows that more of 
the mid-J emission in ULIRGs can be attributed 
to higher bulk excitation of the total molecular gas, 
not just isolated emission from high excitation gas. 

Future work will include the detailed, two-component 
excitation modeling of galaxy spectra with at least eight 
of the thirteen GO lines shown here, as in K14. 
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Fig. A1.— Line Fitting of NGC4388. Top Left: The PSW map of NGC4388 illustrates that this galaxy is extended relative to the largest 
and smallest spectrometer FWHM (white, green circle). As a result, the original spectrum (cyan in Top Right) shows a noticeable gap 
betwee n the SLW and SSW bands. The source/beam size and photometer-matched corrected spectrum is shown in black (described in 
Section [2^. Middle Row, Left Three Columns: Zoomed-in views of the baseline-subtracted spectra (black) for three lines, and the best-fit 
line prome using the FTFitter (red). Bottom Row, Left Three Columns: Probability distribution functions of the integrated line fits of the 
row above, for the FTFitter (red, assuming a Gaussian profile and using the fit and parameter error as median and sigma), and for our new 
method (blue). For some lines, the distribution function is much wider, and/or more tending towards zero than the least squares fitting 
routine would reveal, given the surrounding noise profile. Right Column: Resulting SLEDs for the line fits in Jy km/s {Middle Row) and 
in Lco/^CO,!—0 (Bottom Row). The original fits are shown in red, the SLED we use for fitting from our new method is in blue. 
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Karin Sandstrom for sharing planetary calibration ob- 


APPENDIX 

A. ILLUSTRATED EXAMPLE OF LINE FITTING PROCEDURE 

We chose NGC4388 as an example of a semi-extended galaxy with a fairly good spectrum that degrades in signal/noise 
by the time it reaches the higher-J CO lineS; _ 

As an overview, the top rows of Figure A1 illustrate the source/beam correction described in Section |2.2| The 
photometer PSW map shows that the emission is somewhat extended relative to the SPIRE FTS beam sizes, which 
results in a discontinuity in the original spectrum (cyan, right plot). The corrected spectrum removes this discontinuity, 
and shows the flux emitted in a 43"5 beam. The bottom two rows illustrate high signal/noise (first column) vs. poor 
signal/noise (next two columns) CO lines. Fitting these two lines with a least-squares fitting routine, such as the 
FTFitter, often produce integrated fluxes of signal/noise greater than 3, because the “ringing” in the spectrum is 
well-fit by the intrinsic line profile of the spectrometer. However, inspection of the lines should lead one to question 
why the surrounding ripple peaks are not also high detections of other lines; none of which are expected to be nearly 
the intensity of CO. The resulting probability distribution functions in the bottom rows are thus wider and more 
heavily weighted towards zero. 

To describe how the blue PDFs in the aforementioned figure were created, we focus on the CO J= 12—11 line, for 
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which the FTFitter returns an integrated flux of 0.0298 ± 0.0075 Jy cni“^ centered at 45.73 cni“^. For this procedure, 
we consider the frequency range ± 2 cm from this center, masking out ± one line profile FWHM (0.048 cm“^) 
around any CO, [Cl], and [Nil] lines in this region (in this case, only the CO J= 12—11 line itself). We create a grid 
of injected line amplitudes, /true from 0 to 0.067 (the range is defined by the minimum of 0 or the flux - ba to the flux 
+ 5a). For each amplitude, we inject a sine function of that amplitude at a location within our frequency range and 
refit the spectrum, recording the total measured integrated flux. This procedure is done at evenly sampled frequencies, 
every 0.048 cm“^, over the frequency range (about 80 samples if no other lines are present nearby). For this input 
amplitude, we now have a histogram of measured amplitudes, fobs- All together, we now have a two-dimensional map 
of input fluxes vs. measured fluxes P{fobs\ftrue)j fro m w hich we can back out the probability of the input flux given 
our measured flux. The blue PDF shown in Figure Al is a slice of this map at measured flux of 0.0298 Jy cm“^ 
(P(/trMe|0.0298); in other words, it is the distribution of input fluxes that produced a measured flux of 0.0298 Jy cm“^ 
in this frequency range. 
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